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^ . Abstract 

D ■ In this paper, a Gaussian two-way relay channel, where two source nodes exchange messages with each other 

through a relay, is considered. We assume that all nodes operate in full-duplex mode and there is no direct channel 
between the source nodes. We propose an achievable scheme composed of nested lattice codes for the uplink and 
structured binning for the downlink. We show that the scheme achieves within ^ bit from the cut- set bound for all 
channel parameters and becomes asymptotically optimal as the signal to noise ratios increase. 

Index Terms 

O , Two-way relay channel, wireless networks, network coding, lattice codes 

. . I. Introduction 

OO \ We consider a two-way relay channel (TRC), as shown in Fig. [T] (a). Nodes 1 and 2 want to exchange 
^ messages with each other, and a relay node facilitates the communication between them. This TRC can 
C<l be thought of as a basic building block of general wireless networks, along with the relay channel [[U, 
^ the two-way channel [2J, etc. Recently, there have been a great deal of interest in the capacity of wireless 
O networks. Inspired by network coding [3J, TRC has been studied in the context of network coding for 
^ wireless networks due to its simple structure. However, the capacity region of the general TRC is still 
^ ■ unknown. 

•j-j ' In [4J, several classical relaying strategies for the one-way relay channel [IJ, such as amplify-and- 
rS forward (AF), decode-and-forward (DF), and compress-and-forward (CF), were extended and applied to 
c3 the TRC. AF relaying is a very simple and practical strategy, but due to the noise amplification, it cannot 
be optimal in throughput at low signal to noise ratios (SNRs). DF relaying requires the relay to decode 
all the source messages and, thus, does not suffer from the noise amplification. In [5J, it was shown that 
the achievable rate region of DF relaying can be improved by applying network coding to the decoded 
messages at the relay. This scheme is sometimes optimal in its throughput Q, but it is generally subject 
to the multiplexing loss 

In general, in relay networks, the relay nodes need not reconstruct all the messages, but only need to 
pass sufficient information to the destination nodes to do so. CF or partial DF relaying strategies for the 
TRC, in which the relay does not decode the source messages, were studied in [8]], (91. It was shown 
that these strategies achieve the information theoretic cut- set bound [23 J within a constant number of bits 
when applied to the Gaussian TRC. In [lOJ, [11 J, structured schemes that use lattice codes were proposed 
for the Gaussian TRC, and it was shown that these schemes can achieve the cut-set bound within | bit. 

In this paper, we focus on the Gaussian TRC with full-duplex nodes and no direct communication link 
between the source nodes. Such a Gaussian TRC is shown in Fig. [H (b), and it is essentially the same as 
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(b) 



Fig. 1. Gaussian two-way relay channel 



those considered in [f8]l-[fTT]l. For the uphnk, i.e., the channel from the source nodes to the relay, we propose 
a scheme based on nested lattice codes lfT9ll formed from a lattice chain. This scheme is borrowed from 
the work on the relay networks with interference in fT2], fTT]. By using nested lattice codes for the uplink, 
we can exploit the structural gain of computation coding [15J, which corresponds to a kind of combined 
channel and network coding. For the downlink, i.e., the channel from the relay to the destination nodes, 
we see the channel as a BC with receiver side information (7]], lfT6ll . ifTTll . since the receiver nodes know 
their own transmitted messages. In such a channel, the capacity region can be achieved by the random 
binning of messages [16J. In our strategy, a structural binning of messages, rather than the random one, 
is naturally introduced by the lattice codes used in the uplink. Thus, at each destination node, together 
with the side information on its own message, this binning can be exploited for decoding. 

In fact, as stated above, our work is not the first to apply lattice codes to the Gaussian TRC. However, 
we assume more a general TRC model compared to the other works. In [11 J, it was assumed that the 
channel is symmetric, i.e., all source and relay nodes have the same transmit powers and noise variances. 
In ifTOll . a lattice scheme for the asymmetric Gaussian TRC was proposed. However, the scheme requires 
the existence of a certain class of lattices to achieve a | bit gap to the cut-set bound. This paper extends 
those previous works and shows that we can in fact achieve the cut-set bound within ^ for any channel 
parameters, e.g., transmit powers and noise variances. Moreover, the gap vanishes as the uplink SNRs 
increase. 

This paper is organized as follows. In Section [III we present the channel model and define related 
parameters. The cut-set bound on the capacity region is given in Section |Inl Section JV] illustrates our 
achievable scheme and computes the achievable rate region. Section |V] concludes the paper. 

II. System Model 

We consider a Gaussian two-way relay channel, as shown in Fig. [T](b). We assume that the source and 
relay nodes operate in full-duplex mode and there is no direct path between the two source nodes. The 
variables of the channel are as follows: 

^To be exact, this implies | bit per dimension. For a complex- valued system, as considered in |8|, we have 1 bit gap per complex 
dimension. 
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Vl/^ G {l, . . . , 2^^^}: message of node i, 

T 

: channel input of node i, 

T 

: channel output at the relay, 

T 

: channel input of the relay, 

: channel output at node 

• G {l, . . . , 2^^^}: estimated message of node i, 
where i G {1, 2}, n is the number of channel uses, and Ri denotes the rate of node i. We assume that the 
messages Wi and W2 are independent of each other. Node i transmits ^ at time t to the relay through 
the uplink channel specified by 



(t) 



R 



where Z^^^ is an independent identically distributed (i.i.d.) Gaussian random variable with zero mean and 
variance a^. The transmit signal X-^'' is determined as a function of message Wi and past channel outputs 
y/-^ = |y/^\ . . . i.e., = f^'^ {W^,Y^-^), There are power constraints P„ i G {1,2} on 

the transmitted signals 



t=l 



At the same time, the relay transmits xj^^^ to nodes 1 and 2 through the downlink channel specified by 



y«=xj,nZf^ ZG{1,2} 



where ^ is an i.i.d. Gaussian random variable with zero mean and variance af. The power constraint 
at the relay is given by 



1 ^ 

n ^ 

t=i 



'R 



<Pr^ 



Since the relay has no messages of its own, xj^^^ is formed as a function of past channel outputs 

= {"^R^ • • • .^i'"^^}' i-e-' Xr = Ir (^^"^)- At node 1, the message estimate W2 = Yi) 
is computed from the received signal Yi and its message Wi. The decoding of node 2 is performed 
similarly. Now, the average probability of error is defined as 



For the aforementioned TRC, we say that a rate pair (i?i,i?2) is achievable if a sequence of encoding 
and decoding functions exists such that the error probability vanishes as n tends to infinity. The capacity 
region of the TRC is defined as the convex closure of all achievable rate pairs. 



III. An upper bound for the capacity region 

By the cut-set bound f23], if the rate pair R2) is achievable for a general TRC, a joint probability 
distribution X2, xr) exists such that 

R, < min {/(Xi; Yr, Y2\Xr, X2), /(Xi, X^; FslXs)} , (la) 
R2 < min {/(X2; Yr, Y,\Xr, Xi), /(X2, Xr; Y,\X,)} . (lb) 
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In particular, for the Gaussian TRC, we can use the fact that there is no direct path between nodes 1 
and 2, i.e., (Xi, X2, Yr) Xr (Fi, 1^2), and that Xr (Xi, X2) Yr. This induces another upper 
bound from ([T]), given by 

i?i <min{/(Xi;y^|X2),/(X^;y2)}, (2a) 
i?2 < min{/(X2;yi,|Xi),/(X^;yi)}, (2b) 

for some p{xi, X2, xr). It can be easily seen that, for the Gaussian TRC with transmit power constraints, 
all terms under the minimizations in (O are maximized by the product distribution p{xi,X2,xr) = 
p{xi)p{x2)p{xr), where p{xi), p{x2), and p{xr) are Gaussian probability density functions with zero 
means and variances Pi, P2, and Pr, respectively. The resulting upper bound on the capacity region is 
given by 

ii.<mi„{llog(l + |),llog(l + |)}. (3a) 
i!.<m,„{llog(l + |),ll„g(l + ^)}. (3b) 



IV. An achievable rate region for the Gaussian TRC 

In this section, we compute an achievable rate region for the Gaussian TRC. For the uplink, we consider 
using nested lattice codes, which are formed from a lattice chain. For the downlink, we use a structured 
binning of messages at the relay, which is naturally introduced by the nested lattice codes. The destination 
nodes decode each other's message using this binning and the side information on their own transmitted 
messages. 

The main result of this section is as follows: 

Theorem 1: For a Gaussian TRC, as shown in Fig. [T] (b), we can achieve the following region: 




(4a) 
(4b) 



where = max{x,0}. 

Note that the achievable rate region in Theorem [T] is within | bit of the upper bound dS]), regardless 
of channel parameters such as the transmit powers and noise variances. Moreover, as the uplink SNRs 
^ and p- increase, the gap vanishes and our achievable region asymptotically approaches the capacity 
region of the Gaussian TRC. 

We prove Theorem [T] in the following subsections. 



A. Lattice scheme for the uplink 

For the scheme for the uplink, we consider a lattice coding scheme. We will not cover the full details 
of lattices and lattice codes due to page limitations. For a comprehensive review, we refer the reader to 
[fT9ll - ll2Tll and the references therein. 

A nested lattice code is defined in terms of two n-dimensional lattices Ag< and A^, which form a lattice 
partition Ag</A^, i.e., A^ C Ag.. The nested lattice code is a lattice code which uses Ag< as codewords 
and the Voronoi region of A^ as a shaping region. For Ag</A^, we define the set of coset leaders as 

C = {Ag mod A^} ^ {Ag n 7^}, 
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Fig. 2. Example of a lattice chain and sets of coset leaders. C2 C Ci C Ac- 



where TZ is the Voronoi region of A. Then the coding rate of the nested lattice code is given by 

1 . 1^1 1 . Vol(A^) 
n n Vol(Ag<) 

where Vol(-) denotes the volume of the Voronoi region of a lattice. For the TRC, we should design two 
nested lattice codes, one for each source node. This subsection will show how the nested lattice codes 
are formed. In the following argument, we assume that Pi > P2 without loss of generality. Now, let us 
first consider a theorem that is a key for our code construction. 

Theorem 2: For any A > P2 ^ and 7 > 0, a sequence of n-dimensional lattice chains A^ C A2 ^ Ag. 
exists that satisfies the following properties. 

a) A^ and are simultaneously Rogers-good and Poltyrev-good while Ag< is Poltyrev-good (for the 
notion of goodness of lattices, see [l20l ). 

b) For any e > 0, — e < cr^(A^) < P^, i G {1,2}, for sufficiently large n, where cr^(-) denotes the 
second moment per dimension associated with the Voronoi region of the lattice. 

c) The coding rate of the nested lattice code associated with the lattice partition Ag</A2 can approach any 
value as n tends to infinity, i.e., 

R2 = - log \C,\ = 1 log f J^ij^'l = 7 + 0.(1), (5) 



n ' n ^ VVol(A^), 

where C2 = {Aj^modA2} and o„(l) ^ as n ^ cx). Furthermore, the coding rate of the nested lattice 
code associated with A^/A^f is given by 

where Ci = {A^ mod A^}. 

Proof: See the proof of Theorem 2 in ifTSll . ■ 
For instance, a lattice chain and the corresponding sets of coset leaders are visualized in Fig. [21 for the 
two-dimensional case. 

Encoding 

Let us think of a lattice chain (more precisely, a sequence of lattice chains) and sets of coset leaders 
as described in Theorem [21 We use Ci and C2 for nodes 1 and 2 respectively. For node i, the message 
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set {l, . . . ,2^^^} is one-to-one mapped to C^. Thus, to transmit a message, node i chooses G Ci 
associated with the message and sends 

X, = (W, + U,) modA,, 

where is a random dither vector with \]m.i{lZi) and IZi denotes the Voronoi region of (we 

suppressed the superscript for simphcity). The dither vectors U^, i G {1,2}, are independent of each 
other and also independent of the messages and the noise. We assume that each is known to the source 
nodes and the relay. Note that, due to the crypto-lemma [21 J, is uniformly distributed over TZi and 
independent of W^. Thus, the average transmit power of node i is equal to cr^(A^), which approaches Pi 
as n tends to infinity, and the power constraint is met. 

Decoding 

The received vector at the relay is given by 



where Tar 



7(1) 



R 



Xi + X2 + ^R^ 



. Upon receiving Y^^, the relay computes 



R 



mod Ai 



where 



T = 



2 2 
5](W, + U,)modA,-X;X 

2 2 
a Xj + aT^R — Uj mod Ai 

T + Z^) modAi, 



mod Ai 



^(w,-g,(w, + u,)) 



Li=i 



[Wi + W2 - g2(W2 + U2)] mod Ai, (6) 
-(l-a)(Xi + X2) + aZH, (7) 

a e [0, 1] is a scaling factor, and Qj{-) denotes the nearest neighbor lattice quantizer associated with Aj. 
If we let a be the minimum mean-square error (MMSE) coefficient 

P1 + P2 



a 



the variance of the effective noise © satisfies 



n 



< 



(Pi + A)4 



From the chain relation of lattices in Theorem [2l it follows that T G Ci. Moreover, using the crypto-lemma, 
it is obvious that T is uniformly distributed over Ci and independent of Z^^ [13, Lemma 2]. 

The relay attempts to recover T from Y^^ instead of recovering Wi and W2 separately. Thus, the 
lattice scheme inherits the idea of computation coding [jTSl and physical-layer network coding [181 . Also, 
by not requiring the relay to decode both messages, Wi and W2, we can avoid the multiplexing loss 
at the relay. The method of decoding we consider is Euclidean lattice decoding [fT9l -[l22l. which 
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finds the closest point to Y^^ in A^. Thus, the estimate of T is given by T = Qc {^^^ where Qc{') 
denotes the nearest neighbor lattice quantizer associated with A^. Then, from the lattice symmetry and 
the independence between T and Zr, the probability of decoding error is given by 

Pe = Pr{t7^T} 

= Pr|z^modAi ^7^c|, (8) 



where TZc denotes the Voronoi region of A^. We then have the following theorem. 
Theorem 3: Let 

r Pi , Pi 



Rl = 



log 



R 



For any Ri < Rl and a lattice chain as described in Theorem |2l with Ri approaching Ri, i.e., Ri 
Rl + On{l), the error probabiUty under EucUdean lattice decoding © is bounded by 



Pe < e 



-n(f;p(22(«!-«i))-o„(l)) 



where Ep{-) is the Poltyrev exponent [|22l|. 

Proof: See the proof of Theorem 3 in (13]. ■ 
According to Theorem |3l the error probability vanishes as n ^ oo if < i?]' since Ep{x) > for 
X > 1. This implies that the nested lattice code can have any rate below R^ for the reliable decoding of 
T. Thus, by c) of Theorem |2] and Theorem [3l the error probability at the relay vanishes as n — > oo if 



Ri < 



6 . 6 

P1 + P2 4 



1,2. 



(9) 



B. Downlink phase 

We first generate 2^^^ n-sequences with each element i.i.d. according to A/'(0, Pr). These sequences 
form a codebook Cr. We assume one-to-one correspondence between each t G Ci and a codeword 
^ Cr. To make this correspondence explicit, we use the notation Xi^(t). After the relay decodes T, it 
transmits Xi^(T) at the next block to nodes 1 and 2. We now assume that there is no error in the uplink, 
i.e., T = T. Under this condition, T is uniform over Ci, and, thus, Xi^(T) is also uniformly chosen from 
Cr. 



Upon receiving Yi = X/j + Zi, where Zi 



7(1) 7W 



1 T 



, node 1 estimates the relay message T 



as Ti = ti if a unique codeword exists in Cr^i such that (X/i(ti), Yi) are jointly typical, where 

Cr,i = {XR(t) : t = [Wi + W2 - (32(w2 + U2)] mod Ai, W2 G C2} . 
Then, from the knowledge of Wi and Ti, node 1 estimates the message of node 2 as 

W2 = (±1 - Wi) mod A2. (10) 

2'*^2. Thus, from the 

(11) 

(12) 



Given T = T, we have W2 = W2 if and only if Ti = T. Note that |Cr,i| 
argument of random coding and jointly typical decoding Il23]l . we have 



Pr<^Ti 7^ T|T 







as n 



00 if 



i?2< ^log (1 + ^ 
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Similarly, at node 2, the relay message is estimated to be T2 = t2 by finding a unique codeword in 
Cr^2 such that (XR(t2), Y2) are jointly typical, where 

Cr,2 = {Xfl(t) : t = [wi + W2 - g2(W2 + U2)] mod Ai, wi G . 

Then the message of node 1 is estimated as 



Wi= T2-W2 + g2(W2 + U2) modAi. 



(13) 



Since \Cr^i 



2nRi^ we have 



Pr|t2 7^ t|t = T 








(14) 



as n ^ 00 if 



(15) 



Note that, in the downhnk, although the channel setting is broadcast, nodes 1 and 2 achieve their 
point-to-point channel capacities ([T2h and ([T5]) without being affected by each other. This is because of 
the side information on the transmitted message at each node and the binning of message. In our scheme, 
the relation in ^ represents how the message pair (Wi, W2) is binned to T. 

C. Achievable rate region 

Clearly, the message estimates ([TOh and ([T3]) are exact if and only if Ti = T2 = T. Thus, the error 
probability is given by 



By Theorem [3l the first term of vanishes as n ^ oc if i?^ < i?*, i G {1, 2}. Also, by ([11]) and ([T4h . 
the second and third terms also vanish as n ^ oc if ([T2h and (fTSl) hold. Thus, the achievable rate region 
& follows from ([9]), ([El), and (O. 

V. Conclusion 

In this paper, we considered the Gaussian TRC. An achievable scheme was presented based on nested 
lattice codes for the uplink and structured binning for the downlink. The resulting achievable rate region 
approaches to within \ bit of the cut-set bound for all channel parameters, and the gap eventually vanishes 
in the high SNR regime. Though the capacity region is very nearly reached, the exact capacity region of 
the Gaussian TRC is still an open problem. 
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